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ABSTRACT 

We study the continuous limit of the multiple gravitational lensing theory 
based on the thin lens approximation. Under the approximation, we define a 
new, light-path dependent angular diameter distance d and show that it satisfies 
the optical scalar equation. The distance provides relations between quantities 
used in the gravitational lensing theory (the convergence-, the shear- and the 
twist-term) and those used in the scalar optics theory (the rates of expansion, 
shear and rotation). 

Subject headings: cosmology:theory — distance scale — gravitational lensing 



1. Introduction 



The distance-redshift relation is one of the most important and difficult subjects in 
observational cosmology (here, "distance" means the angular diameter distance). Recent 
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observations on Type la supernovae with high-redshifts suggest the large deacceleration 
parameter q = VLq/2 — A (^o, A are the present values of the density parameter and of the 
cosmological constant, respectively; Perlmutter et al. 1999). The reliability of the suggestion 
depends on that of the distance measure to supernova adopted in the observation. 

Since our universe is assumed to be described by the Friedmann-Lemaire universe which 
is homogeneous and isotropic in a large scale, we have used the Mattig formula as the 
distance measure (e.g. Weinberg 1972, and see also Appendix Al). However we know that 
the universe is not homogeneous in a small scale and that light ray propagating in the locally 
inhomogeneous universe is gravitationally affected by clumps intervening between the source 
and the observer which is known as the gravitational lens effect. Then we have to estimate 
the distance to a source based on the light ray propagating in the homogeneous universe 
in a global scale but locally inhomogeneous universe. Unfortunately, however, we have not 
known the solution of the Einstein field equations describing such a universe yet. 

A recipe to obtain a distance-redshift relation without the exact solution of the field 
equations was presented by Dyer & Roeder (1972, 1973). They gave a distance measure 
based on light bundle propagating away from all the clumps between the source and the 
observer, where the rate of shear in the scalar optics theory is negligible (the shear-free 
assumption). In general, their distance measure is larger than that given by the Mattig 
formula. Consequently, it yields the dimming effect on the observed flux compared with one 
based on the Mattig distance measure. They also showed that the distance measure depends 
on the dimming effect more seriously than the the gravitational magnification effect in the 
Swiss-cheese model (Bertotti 1966; Dyer & Roeder 1974). 

However the shear-free assumption does not always hold for light bundles propagating 
in the general space-time. In fact Watanabe & Sasaki (1990) showed that when the light 
ray experiences multiple gravitational scattering due to intervening clumps, the contribution 
of shear-rate to the distance measure (the Weyl focusing) is comparable with one of matter 
inside the beam (the Ricci focusing). Taking into account the single gravitational lens effect, 
Weinberg (1976) showed that an average flux from source in a clumpy universe with a low 
deacceleration parameter q is equal to the flux in the Friedmann-Lemaire universe with the 
same go- Peacock (1986) developed Weinberg's argument to the case with an arbitrary fl 
and the smoothness parameter a ~ 1. Their results mean that the Dyer- Roeder distance 
multiplied by (fi) ((/i) is the averaged value of the gravitational magnification factor /i) 
can be considered to be equivalent to the distance obtained by the Mattig formula. The 
averaged distance may be useful in analyses such as iV — m or iV — z relations, in which the 
gravitational lensing effect must be statistically taken into account (e.g. Omote & Yoshida 
1990; Yoshida & Omote 1992). 
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However, in the analysis of the distance-redshift relation for individual objects the av- 
eraged distance is no longer available (Hammer & Nottale 1986). Then we are faced with 
the problem how we should estimate the distance to a specific object? We can assume that 
the distance may be given by d = n~ l / 2 d (d denotes the Dyer-Roeder distance) because the 
observed flux from the source is brighter by factor fi than one with no lensing effects in the 
clumpy universe. Although the above definition of the distance d = n~ l t 2 d seems reasonable, 
it is not a trivial problem to show that the distance d to the source is coincident with one 
obtained from the optical scalar equations, which is the subject that we will discuss in this 
paper. We will show that a light-path dependent distance (i~ l ^ 2 d satisfies the Raychaudhuri 
equation for a null geodesic and that the rates of expansion and shear in the scalar optics 
theory can be expressed in terms of quantities of the gravitational lens theory under the thin 
lens approximation. 

In this paper we assume that the universe is described, on average, by the Friedmann- 
Lemaitre model, but locally by the clumpy universe given by Dyer & Roeder (1973). The 
same relation between the cosmological time and the redshift 1 is used in both universe models 
(the homogeneous universe and the clumpy universe). On the other hand, the distances to 
an object from an observer and from another object when gravitational lensing effects are 
absent are assumed to be given, not by the Mattig formula, but by the Dyer-Roeder distance. 

The organization of this paper is as follows. After a brief review of the multiple grav- 
itational lens effect in §§2.1, we give a formulation of the continuous limit of the multiple 
gravitational lens effect (§§2.2). In §§3.1, a light-path dependent angular diameter distance 
d will be defined, and we will obtain equations satisfied by the elements of the Jacobian 
matrix A which gives a mapping from an observer plane to a source plane at redshift z. In 
§§3.2, we will investigate relations between quantities in the gravitational lens theory and 
those in the scalar optics theory. Finally, we shall give the summary and conclusion of this 
paper. 

2. Basic Formulation 

In this section we will present a brief review of multiple lensing effects and give a 
formulation of the continuous limit of the multiple lensing. Furthermore we will derive an 
equation satisfied by the Jacobian matrix of the multiple gravitational lensings. 



1 Usually the redshift is given by k^u^, where and u M are the tangent to a null geodesic and an 
observer's four velocity. Therefore the redshift of a object depends on the distribution of inhomogeneities. 
However, here, we ignore such dependence of the redshift (Dyer & Roeder 1974). 
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We assume that our universe is described by the Friedmann-Lemaire model in the very 
large scale (see Appendix Al), but by the clumpy model (Dyer & Roeder 1972, 1973) in 
the small scale. The fraction of the smoothly distributed matter to the average density p 
of the universe is a (the smoothness parameter), and the rest part of the matter in the 
universe is assumed to be concentrated into the clumps which act as the gravitational lenses. 
In this model the angular diameter distance D{z\\ z 2 ) from a specific object at redshift 
Z\ to another at z 2 is given by the Dyer-Roeder angular diameter distance D(zi,z 2 ) (see 
Appendix A2). We introduce an dimensionless angular diameter distance from Z\ to z 2 as 
d(zi, z 2 ) = D(zi, z 2 )/(c/H ), where H is the Hubble constant. 



2.1. Multiple lens equations 

Suppose that there are N lenses which are randomly distributed at redshift Zi(0 < z\ < 
z 2 < ■ ■ ■ < z N ) in the universe and that z s = z N+ i > z N is a redshift of a source. We set 
lens planes perpendicular to the line of sight at each redshift. The origin on each lens plane 
is located on the point intersected with the line of sight. 

The multi-plane lens equation for the source is given by (see, e.g. Schneider et al. 1992; 
Yoshida, Nakamura & Omote 2004) 

where 6j is the angular position on the i-th lens plane of the light from the source, and ctj 
denotes the deflection angle due to the i-th lens. It is useful to introduce a renormalized 
deflection angle ctj defined by ctj = d(0;Zj)ctj. In equation (1) Oi are recursively expressed 
as follows: 

blN - 1>2 (2) 
9 * = »■ - fo -= 2 < 3 > 

From the above equations we see that Q{ in equation (2) can be regarded as a source at 
redshift Z{ for the foreground lenses. Furthermore, we can reduce a relation between Oi 
and 0i + i as follows: 

— = — = -(1 + Zi)ai(0i), (4) 

Xi Xi+1 Xi—1 Xi 

where Xi is defined in the Appendix A3. The same equation is also derived from the time 
delay of the multiple lensing effects with the Fermat principle (Blandford & Narayan 1986). 
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The deflection angle oli due to the i-th lens in the clumpy universe is given by 

where Azi = z% — Z(z) = c[T — T(z)] (T and T(z) are the present cosmological time 
and a cosmological time when the light emitted from the source at z, respectively) and 
A a (0',Z(z t )) defined by 

8a P (D(0; z)0, Z{z)) = p(D(0; z)0, Z{z)) - ap(z) = p (l + zfA^O, z) (6) 

is an inhomogeneity of the matter distribution on the i-th lens plane. 

From equation (4) it can be shown that the Jacobian matrix A« = d6i/d6i satisfies the 
differential equation 



= -(1 + ^)0,(0,)^, (7) 



Xi Xi+1 X.i—1 Xi 

where Uj is defined as 

Ui(9() = g^!) = 3^ (1 + ^(0;^ ! ft M A ^ z{Zi))(l ' iei _ n (8) 

Ot>i Z Y [Zi) 7T JJ V 

In equation (9) T = (r x , T y ) denotes the shear-term due to the i-th lens, and is given by 

r *^) = r ^) = ^r- ( 10 ) 

Seitz & Schneider (1992) derived a recursive formula of the Jacobian matrix Aj equivalent 
to equation (7). However, their formula is not useful to derive a differential equation of the 
Jacobian matrix by taking a continuous limit, as mentioned below. 



2.2. Continuous limit of the multiple lens effect 



In above equations (4) and (7), we take the limit of Zi+i — > Zi(Azt — > 0) and obtain the 
following differential equations for 0(z) and A(z) 



d_ 

dz 



dz d 
dx dz 



0(z) 



JJD 



d 

dz 



3 (1 + z) 3 d 2 (0;z) 1 

2 F(iT 



^JJj 2 e'A^e',z(z))()'(e(z)-e') 



\d(z) - o' 



A(0(z)). (12) 
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The equation (11) has a formal solution 

r(a o 3 n f ^ (i + C)MQ;CMC;*) i ff 2 0(0-0' fl . 

e(z) = o --n j Q d( — d(M)y(c) nJJ v deA ^°' z{c)) \o(0-er ( 3) 

provided that 0(2) satisfies initial conditions 0(0) = 0o an d d&/dz\ z= o = (see Appendix 
Bl). If the intervening lenses distribution is known, 0{z) can be obtained in the form of a 
power series with Qq. A source at 0(z) in the no lenses case (A g (0', Z(z)) = 0) would be 
detected in the direction to the angular position 0q. 

We can also obtain the following formal solution of equation (12) with the initial con- 
ditions A(0(O)) = 1 (2x2 unit matrix) and dA/dz\ z=0 = O (the 2x2 zero matrix), 

A W ,)) = . - 5* [ ^JM0^oi u{e{O , om , (14) 

where the matrix U is given in terms of its elements as follows : 



(i(0(z) : z) = - II d 2 0'A g (0',z)U'(0(z)-0') 



IT 



V 



A^0(Z),Z)- 1X (0(Z),Z) -ly(0(z),z) 

-l v {0{z),z) A„(0(z),z)+ lx (0(z),z) 



(15) 



A(0(z)) = l-^o [ X(Z) d X i r d X2 (l + C 2 ) 5 rf 4 (0; C 2 )U(0(C 2 ), C 2 )A(0(C 2 ))- (17) 

1 Jy(0) Jx(0) 



In equation (15) ~f(0(z),z) is defined by 

j(0(z), z) = lJJ v d 2 0'A^0>, z)T (0(z) - 0') . (16) 
As described in Appendix B2, equation (14) can be rewritten in the form: 

'x(o) ^x(o) 

Using this expression, we can easily obtain the second derivative of A with respect to x( z ) : 

A"( X ) = -W"( X )A( X ), (18) 

where the matrix W is given by 

q rx{z) rxi 
W( X ) = -fio / d X i ^ 2 (l + C 2 ) 5 rf 4 (O;C 2 )U(0(C 2 )) 

(19) 



r x(o) J X (P) 

<X)-S x (x) Syix) 

-S y (x) k( X ) + S X ( X ) J ' 
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In equation (19) the convergence k(x) and the shear S(x) produced by a single lensing are 
defined as follows: 



q rx(z) rxi 

no/ dXi / rfx 2 (l + C2) 5 rf 4 (0;C 2 )A ci ^(C 2 )), (20) 

q rx{z) rxi 

S{x) = / rfX2(l + C 2 ) 5 rf 4 (0;C 2 )7WC2)). (21) 



In equation (18), A(x), W(x) and the prime ' denote A(0(z(x))), W(0(z(x))) and the 
derivative with respect to Xi respectively. Appendix B2 gives a proof that equation (18) is 
the equivalent to equation (12). In evaluation of A, we must notice that A depends on 0{z) 
which is the solution of equation (11). 



3. Magnification Factor and the Optical Scalar Equations 

3.1. Jacobian matrix 

In the clumpy universe (Dyer & Roeder 1973) if the light ray is not affected by any 
clump which intervenes between a source and the observer, the flux f from the source at 
redshift z with absolute luminosity L would be detected as 

/o = 4nic/H yil + zYd 2 i0;z)- (22) 

However, in the real case, observers detect a light ray which is gravitationally lensed by the 
clumps. Therefore the observed flux / b s is given by / Q b s = fJ>fo, where \i is the gravitational 
magnification factor defined by fx — det A -1 . It is important to notice that \x depends on the 
light-path. Since observers recognize that the angular diameter distance to the source with 
flux / obs is given by (H /c) x L 1 / 2 [47r(l + -2) 4 / bs]~ 1 ^ 2 ) the new angular diameter distance d 
is defined by 

d(0(2)|O; z) = /2- l/2 i0iz))diO; z). (23) 

We should notice that d has a dependence on the light-path. Then, even if two sources have 
a same redshift and a same absolute luminosity, their distances are not always reduced to 
be same values. In the below we investigate the equation of the new distance d. Hereafter 
we assume /i ^ 0, i.e. the light ray considered here does not pass through a conjugate point 
(Wald 1984). 

For this purpose we need, first, to obtain an equation of the gravitational magnification 
factor /i. In general, Jacobian matrix A from the observe plane to a source plane at redshift 
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Kx + G X Gy + Ky 



Go, 



Ky K x 



z is expressed as 
A = 

where oVs are the Pauli matrices: 



G x 



K x a + G^ai + iK y a 2 + £^(73, 



(24) 











1) 


, oi = ( 



-i 

1 



, &3 = 



1 

-1 



(25) 



In this expression, the magnification factor ji is given by n = det A" 1 = (K K - G 
G)^ 1 , where K x ,K y and G denote the cumulative convergence-, twist- and shear-terms, 
respectively. In a single lensing case, A is symmetric, i.e., the twist-term K y always vanishes. 
On the other hand, in the multiple lensing case, it does not vanish, in general. 

According to Seitz et al. (1994), we can define an optical deformation matrix Q in terms 
of A and its derivative with respect to x as follows: 



which gives 



Q = A' A" 1 , 



QA. 



(26) 
(27) 



Inserting equation (24) into equation (26), we obtain an explicit form of the optical defor- 
mation matrix: 



Q 



+ie ab (k a k' b + G' a G b ) a 2 + 5 ab (k a G' b - K' a G b ) a 3 , 



(28) 



where e xx = e yy = 0, e xy = -e yx = 1, K a = fi 1/2 K a and G a = ^ 1/2 G a . 



Using a relation (18), the definition of Q and its derivative with respect to x, we can 
obtain a relation between Q and W" as follows 



Q' + Q 2 = -W", 



(29) 



where W" is written in terms of the Pauli matrices as W" = k,"<j — SyG\ — S"cr 3 . Using 
equations (28) and (29), we find equations of the coefficients of the Pauli matrices in Q in 
the following 



(ln / ,-V2)" + (l n/i 



-l/2\' 



e ab (k a G' b - k' a G b ^j 
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+ 

Sab 



5 nh ( K„G' h - K'G, 



«i> 1 K a G' h - K' a G b 



£ab 



i, ~ K' a G h 



e ab [ K a G' b - K' a G b 



e ab [K a K' b + G' a G 

+ 2{\n^)'[5 ab [k a G' b 
+ 2(ln^ 1 /2)' r 

(k a K' b + G a G' b )} ' + 2 (ln/i" 1/2 )' [e ab (k a K' b + G a G' b ) 



on 

Qll 

>=>y, 
0. 



(30) 

(31) 
(32) 
(33) 



Equation (30) is equivalent to one derived by Seitz et al. (1994). They gave the equation 
of the magnification factor in the clumpy universe by using a rate of shear derived from a 
Weyl term of the scalar optics theory in the linearized general relativity. On the other hand, 
equation (30) is derived by taking the continuous limit of the multiple gravitational lens 
theory. It should be noticed, therefore, that in our formulation we can obtain the right hand 
side of equation (30) from data of lens distributions itself, without referring to the Weyl term 
(the metric terms given by solving the Einstein equation). 

The second derivative of d with respect to the afline parameter v can be expressed as 
combination of the second derivative of d with respect to v and the second derivative of p^ 1 ! 2 
with respect to x- 

,f #-" I (»W). (34) 



Combining equation (30) with equation (34), the Dyer-Roeder equation (A13) yields 



di 
where 



|-e + e 2 + (s ab t ab y+ ^ a6 s a6 )'-(£ a 5^5) 2 = -^p( J D(o^)0^),z(^))(i + ^) 2 , (35) 



e 

&ab 



d 

dv 



H 



\nd(6(z)\0;z), 



K a G' b - K' a Gb 
d 2 (0;z) ' 

k a k' b + G a G' h 



d 2 (0;z) ' 

Furthermore, inserting equations (37), (38) into equations (31)-(33), we find that 

— (^5 ab Tj ab j + 2Q5 ab Y, ab 

' ' ] + 2Qe ab t a 
dv \ / 

and 



H 2 

d ( ~ \ ~ ~ AtxG 
^ \ £ab^ab J + 2<ae ab h ab = -jjY 



-p(z)(l + z) 2 lx (0(z)), 
p(z)(l + z) 2 ly (0(z)), 



(36) 
(37) 

(38) 

(39) 
(40) 
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d_ 

dv 



(sabtiab) + 2Qe ab uj = 0. 



(41) 



Here we define £ = (S a b + i£ a b)^ab and uj = e a bU a b and then rewrite equations (35), (39)-(41) 
as follows: 



where 



n 



AttG 



H 2 



^-Q + Q 2 + \t\ 2 -Co 2 = Tl, 
dv 

d ~ 

— £ + 20S = T, 

dv 

^—Cj + 2Blj = 0, 
dv 



p(D(0;z)e(z)\0;z)(l + z) 2 , 



(42) 
(43) 
(44) 

(45) 



AnG 

AG 
Hi 



p(z) [ lx {0{z))+i ly {0{z))]{l + zf 



I J J d 2 6'5 a p(D(0; z)0', Z{z))T{0{z) - 0>)(1 + z) 2 , (46) 



and F(z) = V x (z) + iT y (z) denotes the complex shear term due to a lens at redshift z. 

From equations (36) and (42), we find that our new distance (23) satisfies the following 
equation: 



2 4ttG 



HI 



p(D(0;z)e(z),Z(z))(l + z) 



d(0{z)\0;z)=0. (47) 



We should notice that we cannot obtain solution of equation (47) without solving equation 
(11). 



3.2. Relations between the Jacobian matrix and the optical scalars 

Since the cross-sectional area of the light bundle, A, is proportional to square of the 
angular diameter distance (Dyer & Roeder 1972), we can rewrite 6 as 

e = IA lnA (48) 

which is the rate of expansion in the scalar optics theory. Furthermore, equation (45) can 
be regarded as the Ricci term, which is given by TZ = — (1/2) R^k^k". Therefore, we can 
show that equation (35) or (42) can be regarded as the optical scalar equation of the rate 
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of expansion, if S and Q correspond to the rates of shear and rotation in the scalar optics 
theory. 

First, let us investigate Co. Combining initial conditions of A and definitions of 9 and u, 
we obtain a solution of (44), u> — 0. In the scalar optics, the rotation rate oo = ^Juo a poj al3 /2 
always vanishes, because the covariant derivative of the tangent vector of the null geodesic is 
symmetric, i.e., uj a p = V[ a k^ = 0. From this fact and equation (44), we can regard uj as the 
rotation rate in the scalar optics theory. Furthermore, we find that the optical deformation 
matrix Q is symmetric. In a discrete multiple gravitational lensing theory, the symmetry of 
a matrix corresponding to the optical deformation matrix was shown by Seitz & Schneider 
(1992) by means of their recursive formula of the Jacobian matrix. However, in their paper, 
the reasons of the symmetry was not clarified. 

Next, let us notice the right hand side of equation (43), T . This is similar to the Weyl 
term in the scalar optics theory. The difference between T and the Weyl term in the linear 
perturbation theory of the general relativity (Sasaki 1987; Futamase & Sasaki 1989) exists 
in the following point: while the Weyl term at arbitrary point (-D(0; z)0(z), Z(z)) on the 
light-path in the linear perturbation theory has some contribution from the matter in the 
whole universe outside the light bundle, T in our formalism has a contribution from only 
the matter on the lens plane with the same redshift z. This comes from the fact that our 
formalism is based on the thin lens approximation. Therefore, we can regard JF as the Weyl 
term in scalar optics theory as long as the thin lens approximation is valid. 

Finally we find that equations (42)-(44) are equivalent to the optical scalar equations 
(All) in an inhomogeneous universe under the thin lens approximation. Consequently the 
optical scalars are expressed in terms of the elements of the Jacobian matrix and the Dyer- 
Roeder angular diameter distance as follows: 

9 = <L In [fi-^diO; z)],E = - (5 ab + ie ab ) 
where the last equation is exactly zero. 

4. Summary and Conclusions 

In this paper we have formulated the continuous limit of the multiple gravitational lens 
effect by taking the limit in which the number of lens planes is infinite and then the redshift- 
intervals of lenses are infinitesimal with keeping the present average density po fixed. It, 
however, does not mean that the number of lenses must be infinite. In this case, we can 



K a G' b 



K'G a 



d 2 (0;z) 



V = -tab 



k a k' b + G' a G a 



d 2 (0;z) 



(49) 
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interpret that the finite number of lenses have some finite masses in their lens planes and 
other infinite number of lenses have zero masses in the planes, and that in the zero-mass 
lens planes, the deflection angles, d, should be and the matrices, U, should be I. The 
formalism is able to be adopted in the case where there are some finite number of lenses 
between sources and us. 

Using the formalism we found that the new angular diameter distance, d = p~ l l 2 d sat- 
isfies the optical scalar equation of the expansion rate in an inhomogeneous universe. It 
is interesting that the influence of inhomogeneities on the angular diameter distance is com- 
pletely recovered by multiplication of p~ l l 2 to the Dyer-Roeder distance, which is obtained 
by ignoring the influence of all clumps on the light beams. 

We also found that quantities defined in the scalar optics theory (the rates of expansion 
G, shear E and rotation a;) are expressed in terms of quantities defined in the gravitational 
lens theory (the terms of convergence K x , twist K y and shear G). Seitz & Schneider (1994) 
showed that the absolute square of the shear rate, |S| 2 , can be expressed by some combination 
of the determinants and traces of the Jacobian matrices in the discrete multiple gravitational 
lensing theory. It is equivalent to one in our formalism, which is also written as 



This expression gives an immediate estimation of the shear rate by tracing the lens equation 
(11) from the raw data of local distribution of matters. 

The formalism includes the framework of the weak gravitational lensing theory (e.g. 
Bartelmann & Schneider 2001), in which the smoothness parameter a is assumed to be 
unity(the distance measure is given by the Mattig formula and only single gravitational 
lensing is taken in account). In their framework, when the light ray propagates through some 
underdense region, p < p, the light is assumed to be deflected by a lens with some negative 
mass (deflecting mass bp = p — p < 0). Therefore such lenses with negative masses give some 
brightening effect similar to the usual clumps on the flux when the light beam passes away 
from such lenses, because the shear rate due to the negative mass lenses contributes to the 
Raychaudhuri equation in the form of |<5p| 2 . Then there is some possibility that the weak 
lensing theory does overestimate the distance measure due to the negative lenses unless the 
convergence due to the negative mass lens in the light beam balances with the shear outside 
the light beam due to the other ones. In our formalism , unlike the weak lensing theory, 
we do not need to introduce such lenses with the negative masses in order to take into 
account the underdense regions in which the matter density is only considered to be p = ap. 
Therefore there is not the above possibility because the shear in such underdense region does 
not contribute to the distance measure. 



£ = 



2 



G G - K K 

d 4 (0;z) 



Finally, we comment on equation (47) which is the generalized Dyer-Roeder equation. 
Schneider et al. (1992) also mentioned a similar equation, which is slightly different from 
ours. In their equation, the density of clumps b&p = p — ap is not considered, because they 
assumed that the light beam does not pass through bound clumps. However, among real 
situations, there are cases in which the contributions of both the matter density of the bound 
clump in the light bundle and the shear rate of the light bundle are compatible (Watanabe 
& Sasaki 1990). Even in such intermediate case, equation (47) gives the angular diameter 
distance to a target source, by solving the lens equations (11) based on observed distribution 
of foreground clumps. 
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A. Relations between redshift z, afRne parameter v and x-function 

In this appendix, we briefly describe the relation between z (redshift), v (afline parame- 
ter) and x (x - f unc tion). 



A.l. Background universe 



We assume that the universe is described by the Friedmann-Lemaitre universe on aver- 
age. Then the geometry on average is expressed as 

dr 2 



ds z = c 2 dT z - a\T) 



1 — kr 2 



+ r 2 (de 2 



sin 



(Al) 



where T is the cosmological time, a(T) denotes the scale factor of the universe which has the 
dimension of distance and k is related with the curvature of the universe. In this geometry, 
the Einstein field equations can be written as 



Ac 2 kc 2 
? 
kc 2 



8nG_ 
3 3 or 



(A2) 
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where p and A are the mean density of the universe and the cosmo logical constant, respec- 
tively, and the dot denotes the derivative with T. Since a(T) and p are given in terms of their 
present values a and p by a(T) = a /(l + z) and p = p (l + z) 3 (in the matter-dominated 
era), respectively, equations (A4) and (A3) are written as 



a 
a 

da 
d 2 



El [fi (l + zf + A - K(l + zf] = HlY 2 ' 
Q (l + zf/2 - A 



q(z) = 



Y 2 (z) 



(A4) 
(A5) 



where H is the present Hubble constant, f2 = 8ttGp /3Hq, A = Ac 2 /3Hq and K = 
k(? ja\H.Q. The above equations gives relations between f2 , A and K: 



K = n + A - 1, 
g = — Aq, 



(A6) 
(A7) 



where g is the present deacceleration parameter defined as — (aa/d 2 ) . It follows from 
equation (A4) that the relation between T and z is given by: 



-cdT = 



dz 



H (l + z)Y(z) 



(A8) 



Furthermore, in the Friedmann-Lemaitre universe, an affine parameter v defined along 
a null geodesic is given in terms of redshift z by 



dv = 



dz 



(\ + z yY(zy 

In this universe, the angular diameter distance D FL (z 2 ; z±) from z\ to z 2 is given as 



(A9) 



c 1 

D FL {z 2 ; Zl ) = a(T 1 )r{z 2 ; z x ) = — sin 

-no VK(1 + Zi) 



dz 
Yjz) 



(A10) 



This is the Mattig formula. 



A. 2. The Dyer-Roeder distance 

Sachs (1961) gave geometrical equations of light bundle propagating in the general 
space-time from the geodesic deviation for the light ray as follows: 

^-9 + 2 + ISI 2 = K, ^-£ + 2£0 = .7 r , (All) 

dv dv 
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where 6 and £ denote the rates of expansion and shear factor of the light bundle, 1Z and T 
are the Ricci and Weyl terms. The rate of expansion G is defined as = rfln \fA/dv, where 
A denotes the cross-sectional area of the light bundle. Since the angular diameter distance 
D to a source at redshift z is proportional to vC4, equations (All) are rewritten as 

-PI') A ^ (A12) 

In general it follows that the angular diameter distance depends not only on the Ricci term 
but also on the Weyl term, where the Weyl term is due to the inhomogeneities outside the 
light bundle while the Ricci term originates from the matter density in the light bundle, 

Dyer & Roeder (1972, 1973) assumed that for the light beams far away from all the 
inhomogeneities, the contribution of E to the distance D from the source is negligible and 
that the matter-density in the light bundle is given by ap, where a (smoothness parameter) 
is the fraction of the density smoothly distributed in the universe to the mean density of the 
universe, and derived the angular diameter distance in the clumpy universe, which satisfies 
the following equation: 

^D+*£fap(l + z) 2 D = 0. (A13) 

The relation (A9) between the afline parameter along the light ray and redshift allows us to 
express equation (A13) as follows: 

d(z ; z) + jfooc ^ + f d{z ; z) = 0, (A14) 
2 Y{z) 



d_ 

dz 



1 + 



where d(zo] z) denotes the dimensionless angular diameter distance D(zo] z)/(c/H ) from an 
object at redshift z to another at z(z < z ). The initial conditions at z = z are given by 



d(z ; z 

d 

dz 



= 0, 

2 = 20 

*"l, = (i + „V W ' (A15) 



(Dyer & Roeder 1973), where the last condition is the Hubble law at redshift z . In the case 
with a — 1, equation (A10) is a solution of the Dyer- Roeder equation. 



A. 3. The x-function 

The x-function is introduced by Schneider et al. (1992) as follows 

f°° dz 

* Ws / (TTzm^mry (A16) 
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where d(0; z) is the dimensionless Dyer-Roeder angular diameter distance. The relation 
between x, v and z is then expressed as 



dx 



dv 



dz 



d 2 (0;z) (l + z) 2 Y(z)d 2 (0;z)' 



(A17) 



Accordingly we can use the x-function as a variable instead of redshift z or afhne parameter 
v. 

The deference between values of x & t redshifts Zi and at Zj < Zj) has an interesting 
relation to the Dyer-Roeder angular diameter distance: 



(A18) 



In the rest part of this subsection we will prove this relation. First of all, we investigate 
the differential equation which d(0; z)x{ z %i z ) satisfies. Let f(z) denote d(0; z)x{ z i\ z) and 
differentiate it with respect to affine parameter v: 



d 
dv 



d 

dv 



d(0;z) 



X{zu z) + ■ 



dz/dv 



d 
dv 



d(0;z) 



{l + z) 2 Y(z)d(0;z) 
The second derivative of f(z) with respect to v is also given by 

4ttG 



X(zi]z) + 



d(0;z) 



(A19) 



cP 
dv 



2 m 



dv 



2 d(0;z) 



X(zi] z) 



HI 



ap(z)(l + zYf(z). 



(A20) 



Therefore we find that f(z) is a solution of the Dyer-Roeder equation (A13). 



Next we investigate what initial conditions f(z) satisfies. It is evident from equations 
(A18) and (A19) that /(^) = and 



dz/dv\ z=z 



(A21) 



h {l + ZifYiz^-Zi) d(0; Zi y 

From equations (A20) and (A21), we find that (1 + Zi)d(0; Zi)f(z) is equivalent to d(zi, z). 

The proof of equation (A18) is thus completed. Functions similar to the x-f unc ti° n 
are defined in different forms by several authors (Nottale & Hammer 1984; Padmanabhan 
& Subramanian 1988) for their own purposes. Nottale & Hammer defined a distance as 
an optical distance, which is equivalent to the inverse of the x-f unc tion. Padmanabhan & 
Subramanian defined the function as an independent solution of the usual Mattig distance 
which satisfies the Raychaudhuri equation of the null geodesic in the homogeneous universe 
model. 



-17- 



Some Proofs 



B.l. Equation (13) is a solution of equation (11) 

In this appendix, we prove that in the continuous limit, the angular position 0(z) of 
images in the lensing plane at redshift z satisfies equation (11). First of all, we differentiate 
equation (13) with respect to redshift z and then obtain 



(Bl) 

After multiplying dz/dx(= —(1 + z) 2 d 2 (0; z)Y(z)) to both sides of this equation, we differ- 
entiate it with respect to z, again. 



d 

dz 



dz d 
dx dz 



0(z) 



-Mi + zy<r(p;z 



d (d(C;z) 



3 o 



dz \ d(0; z) 

l - J J d 2 e'^{e\z{z)) 
i + C)M0;C)^ 



x- 



0{z) - 0' 
\0(z)-0'\ 2 



no 



dz 



l + z) 2 d 2 (0;z)Y(z 



d (d(C;z) 



x 



^JJj 2 o'A a (e',z(c))- 



dz \d(0;z) 

e(0 - & 



\0(O-er (B2) 

Here, we apply the initial conditions of the Dyer-Roeder angular diameter distance, (A15), 
to the first term of the right hand side of equation (B2) to obtain 



d (d(C;z) 



dz \ d(0; z) 



(l + z)Y(z)d(0;z) 



(B3) 



In the second term of the right hand side of equation (B2), equation (A9) provides that 



d(C;z) 



1 d \ 2 d 

(l + z) 2 Y(z)fo\ ^ 



d(0;z)_ 
d((;z 



d(0;z) 



d(0;z)^d((;z)-d((;z)^d(0-,z) 



(B4) 



(l + z) 2 Y(z) 

vanishes because both d(0; z) and d((; z) are solutions of the same equation (A13). Finally 
we find that equation (B2) is rewritten as 



d 

dz 



dz d 
d\ dz 



0{z) 



3 (l + z) 3 <P(0 ;z)l 
2 Y[zT 



±-JJj 2 d'A^0',Z(z)) 



0(z) - 0' 
0{z)-e'\ 2 ' 



(B5) 



Similarly we can prove that equation (14) satisfies equation (12). 
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B.2. About equation (17) 

In this appendix, we prove that if arbitrary function g(z) has a finite value at z — 0, 

dt il + Z d{ ^ C,Z) 9(0 = r iZ) d X i f X{Cl) d X2 (l + C 2 ) V(0; ( 2 )g(( 2 ), (B6) 



d(0;z)Y(z) ■ ' J x(0) J x{0) 

where xi,2 = x(Ci,2) defined in the Appendix A3. This relation is used in the text when 
equation (17) is derived from equation (14). 

To begin with, we replace d((; z)in the left hand side of (B6) by the ^-function by using 
relation (A18): 

C d(0;z)Y(C) 9{Q 



J o dai+(Td' 2 (o;c r^ Y{ ^~ j g(o 
I dcm fl±g^^o - *m £ .c (1+ g;; (0;C) ,(c) 

»x(z) rx{z) 



+ Ci) 5 ^ 4 (0; Ci)^(Ci) + X(z) / + Ci) 5 ^ 4 (0; GMG). (B7) 

x(o) ^x(o) 



In equation (B7), we introduce a function 

h(x)= [ X rfxi(i + Ci) 5 ^ 4 (0;Cik(Ci), 



and rewrite equation (B7) as 



-^ (i + C) 2 d(Q;CMC;*) ^ w v [ x(z) , d 

= X(0)MX(0))+ T dXiMxi), (B8) 
^x(o) 

where we define x(0)Mx(0)) as a value of x(e)Mx( e )) in the limit of e — > 0. The behaviors 
of x( e ) an d Mx( e )) for small e are x( e ) ~ 0(1/ e) and /i(%(e)) ~ e 3 g(e). We thus find that 
x(e)/i(x(e)) ~ e 2 g(e) — > if cjr(0) is finite. Finally we have completed the proof of equation 
(B6). 
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